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The highest precision in direct mass measurements is obtained with Penning trap mass spectrom- 
etry. Most experiments use the interconversion of the magnetron and cyclotron motional modes 
of the stored ion due to excitation by external radiofrequency-quadrupole fields. In this work a 
new excitation scheme, Ramsey's method of time-separated oscillatory fields, has been successfully 
tested. It has been shown to reduce significantly the uncertainty in the determination of the cy- 
clotron frequency and thus of the ion mass of interest. The theoretical description of the ion motion 
excited with Ramsey's method in a Penning trap and subsequently the calculation of the resonance 
line shapes for different excitation times, pulse structures, and detunings of the quadrupole field has 
been carried out in a quantum mechanical framework and is discussed in detail in the preceding 
article in this journal by M. Kretzschmar. Here, the new excitation technique has been applied with 
the ISOLTRAP mass spectrometer at ISOLDE/CERN for mass measurements on stable as well as 
short-lived nuclides. The experimental resonances are in agreement with the theoretical predictions 
and a precision gain close to a factor of four was achieved compared to the use of the conventional 
excitation technique. 
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I. INTRODUCTION 

The mass and its inherent connection with the atomic 
and nuclear binding energy is an important property of 
a nuclide. Thus, precise mass measurements are eminent 
for various applications in many fields of physics [l|, Q • 
The required precision of the atomic mass depends on 
the physics being investigated. For radionuclides, which 
often have half-lives considerably less than a second, it 
ranges from 5m /m = 10~ 5 to below 10 -8 and for stable 
nuclides even down to Sm/m — 10 -11 . Since the interest 
in masses of short-lived and stable nuclides arises from 
a wide range of applications with different requirements 
on the accuracy, a continuous development of e.g. new 
ion detection and preparation techniques is carried on at 
several facilities world-wide @ . 

In 1949 N.F. Ramsey improved the molecular-beam 
magnetic- resonance method of I.I. Rabi by applying 
oscillatory fields to the transvcrsing molecular beam in 
spatially separated regions @ . This replacement of the 
uniformly applied field led to a linewidth reduction to 
60% and thus to a more precise determination of the 
resonance frequency. For a detailed description see [f| 
and references therein. 

In high-precision mass spectrometry the idea to use 
time-separated oscillatory fields to manipulate the radial 



* Correspondence to: gcorge@uni-mainz.de 



motions of confined ions in a Penning trap was first put 
forward at ISOLTRAP by G. Bollen and coworkers in 
1992 |1 and later tested at the SMILETRAP experi- 
ment [l(| • At this time the correct theoretical description 
of the observed line-shapes was not available. Instead, 
the observations were discussed qualitatively in terms 
of the Fourier transform of the applied pulse sequence. 
However, a study of the precision gain as well as on- 
line mass measurements with the new excitation scheme 
were not performed. In the preceding paper [Tl| a the- 
oretical description of the obtained lineshape has been 
developed in a quantum mechanical framework. In the 
following, the corresponding experimental investigations 
and the gain in linewidth and precision of the "Ramsey 
method" will be presented. 



II. THEORETICAL OVERVIEW 



For a quick orientation we review here in a non- 
technical manner the physical assumptions forming the 
basis of the theoretical model, the approximations neces- 
sary for analytical calculations, and the results that are 
most important for mass spectrometry. Details are found 
in the preceding article [llj . 



2 



A. Penning traps with azimuthal quadrupole 
excitation 

1. The ideal Penning trap 

The electromagnetic field configuration of an ideal Pen- 
ning trap consists of a strong homogeneous magnetic field 
Bo = Boe z in the axial direction and an electrostatic 
quadrupole field E = — V<f>o derived from the potential 
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y ). The motion of a 



single ion of mass m and electric charge q in this elec- 
tromagnetic field is determined by the field parameters 
Bo and U , or equivalently, by the cyclotron frequency 
v c = ll> c /2tt = g£>o/(27T77i) and by the axial frequency 
v z = lo z /2n with lo z = \J AqU /(2zq + Tq), where ro and 
zq denote the inner radius of the ring electrode and the 
distance of the endcap electrodes from the trap center, 
respectively. The dynamics of a single ion in an ideal 
Penning trap is described by three uncoupled harmonic 
oscillators: the oscillator of th e cyclotro n motion with 
the frequency uj + = \(lo c + \J^1 — 2wf ), the inverted 
oscillator of the magnetron motion with the frequency 
LO- = \{u c — — 2to z ), and the axial oscillator with 
the frequency tu z . It is u seful to introduce the abbre- 
viation lo\ = \JuJc — 2lo z , so that io± = \{u c ± 
With appropriately chosen canonical coordinates qk, Pk 
(k = +, — , 3) the Hamiltonian for the motion of a single 
spinless ion can be written [ll| as 
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Canonical coordinates offer an easy access to the quan- 
tized version of the theory. Annihilation and creation 
operators for the oscillator quanta are defined by 

flfe = -7= (qk + ipk) , 



(q k -iPk); (k = +,-,3) (2) 



with commutation relations [a±,aj_] = 1, [a±,a^] = 0, 
and [a±, a^] = 0. The Hamiltonian then takes the form 



Ho = faLa+(a} + a+ + i) — ?kj_(a^a_ + i) 
+tkv z (ala 3 + \). 



(3) 



The quantized version of the theory provides a clear pic- 
ture of the energy level scheme associated with the ion 
motion in an ideal Penning trap (see Fig. [1]). Therefore, 
it is valuable for the identification of the interaction with 
the external rf-quadrupole field that is used to convert 
the magnetron motion into the cyclotron motion. 

Real Penning traps generally possess some anhar- 
monic perturbations that introduce nonlinear terms into 
the equations of motion and couple the three motional 
modes, thus preventing us from finding analytical solu- 
tions for the ion motion. This problem can be minimized 
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FIG. 1: Energy level scheme of the harmonic oscillators for 
a spin-less charged particle in an ideal Penning trap. lu+ is 
the modified cyclotron angular frequency, lu z is the angular 
frequency, and u_ is the magnetron angular frequency. n+, 
n z , and n- denote the corresponding quantum numbers. The 
total energy is given by the sum of the energies of the three 
independent harmonic oscillators. The contribution of the 
inverted magnetron oscillator is negative. Zero point energies 
of the oscillators have been subtracted. 




FIG. 2: Radial segmentation of the ring electrode (top view) 
for the application of radiofrequency fields. A predominantly 
quadrupolar field can be generated by applying a radiofre- 
quency between pairs of opposing electrodes of the four-fold 
segmented ring electrode. 



by avoiding large distances of the ion from the trap cen- 
ter. The theoretical discussion assumes also that all an- 
harmonic perturbations and couplings to the axial oscil- 
lator mode can be neglected. Therefore the description 
concentrates exclusively on the cyclotron and magnetron 
motional modes. 
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2. The ideal Penning trap with azimuthal quadrupole 
excitation 

The application of a radiofrequency potential to the 
four-fold segmented ring electrode (see Fig. \2§ adds to 
the Hamiltonian Q new terms that are periodic with the 
driving frequency Vd — uJd/^TT and have a phase \d- The 
leading term in a multipolar expansion is the quadrupole 
contribution oc (x 2 — y 2 ) cos(u>dt + Xd)i higher multipolcs 
must be minimized as they introduce non-linearities into 
the equations of motion. In the quantized version of the 
theory x and y are components of the position operator 
of the ion and can be expressed in terms of the creation 
and annihilation operators of the modified cyclotron and 
magnetron motional modes. It then becomes apparent 
that the quadrupole contribution actually describes three 
physical processes: Most importantly the absorption of 
an electromagnetic field quantum hujd with simultane- 
ous conversion of a magnetron oscillator quantum Tilo- 
into a quantum of the cyclotron oscillator hw+ , together 
with the reverse transition. This process is dominant if 
the driving frequency v& approximately satisfies the reso- 
nance condition huj^ ~ thuj+ + — hio c . It is described 
by an additional term in the Hamiltonian 

H^t) = hgL- i{ " dt+Xd) aX{t)a-(t) 

+e +i ^ dt+Xd) ai(t)a + (t)^ . (4) 

The real coupling constant g has the physical dimension 
of a frequency It is proportional to the amplitude of the 
rf-quadrupole field, but its value is also influenced by de- 
tails of the trap geometry. As shown later, it determines 
the conversion time r c = ir/2g, which is defined as the 
time required for the full conversion of a state of pure 
magnetron motion into a state of pure cyclotron motion 
by a quadrupole field at the resonance frequency ui c (see 
final results in Eq. ([9|)). For each magnetron quantum 
huj- that is annihilated (created) a cyclotron quantum 
huj + is simultaneously created (annihilated). The inter- 
action obviously conserves the total number of quanta 
iVtot = N + + iV_ present in the system and has no sim- 
ple description in a purely classical picture. The other 
two physical processes described by the quadrupole con- 
tribution are the absorption (emission) of a field quan- 
tum hujd with simultaneous creation (annihilation) of two 
quanta of the modified cyclotron oscillator (driving fre- 
quency Ud ~ 2uj + ), and the absorption (emission) of a 
field quantum hujd with simultaneous annihilation (cre- 
ation) of two quanta of the magnetron oscillator (driving 
frequency uid ~ 2cj_) [l2j]. These two processes are neg- 
ligible if Ud ~ to c . 

The ideal Penning trap with quadrupole excitation is 
now described by the total Hamiltonian H(t) = Hq + 
H\{t). The resulting Heisenberg equations of motion 
for the operators a+(t) and a_(i) are linear and time- 
dependent, they permit a general, exact solution. The 
equations and their solutions correspond closely to those 



familiar from the study of magnetic-resonance or the 
quantum theory of two-level systems. An important re- 
sult is that the interconversion of the magnetron and 
cyclotron motional modes by a quadrupole field of fre- 
quency Vd is periodic with the 'Rabi frequency' Ur = 
w r /{2-k) = ^(g/n) 2 + {v d - v c ) 2 . 

3. Ion trajectories 

While the viewpoint of quantum mechanics was very 
helpful for the identification of the relevant interaction, 
it is in general sufficient for the application to mass spec- 
trometry to consider the ion motion as a classical motion, 
i.e. following trajectories of macroscopic scale. To work 
out this aspect expectation values of the quantum me- 
chanical operators are taken with respect to quasiclas- 
sical coherent oscillator states. Thus, the annihilation 
operators a+(t) and a_(i), obtained as solutions of the 
Heisenberg equations of motion, are translated into two 
complex functions a+(t) and ct-{t) that are denoted as 
the 'complex oscillator amplitudes' of the cyclotron and 
magnetron oscillators at time t. The complex conjugate 
functions a+(t) and a*_{t) correspond to the respective 
creation operators. The explicite solution of the initial 
value problem is [Tlj |: 
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(5) 



(6) 



In these equations 6 — u>d — u c is the detuning of the 
driving quadrupole field, lu r — \J Ag 2 + S 2 the Rabi fre- 
quency, Xd the phase of the quadrupole field at time 
t = 0, and a±(0) = |a±(0)| exp[=Fix±] the initial values 
of the complex oscillator amplitudes. For example, for an 
initial state of pure magnetron motion or pure cyclotron 
motion they are <x+{0) = or a_(0) = 0, respectively. 

The practical importance of the complex oscillator am- 
plitudes for our applications lies in their relation to the 
instantaneous radii for the cyclotron and the magnetron 
motion, 
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(7) 

Relating the complex oscillator amplitudes to the original 
Cartesian coordinates and velocities, an explicite para- 
metric representation of the ion trajectories in the xy- 
plane can be derived 
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A graphical representation of an ion trajectory calcu- 
lated by this approach is shown in Fig. [3J For a three- 
dimensional representation one has to add, of course, the 
oscillatory axial motion. 
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FIG. 3: Conversion of pure magnetron motion to pure cy- 
clotron motion by the action of an azimuthal quadrupolar 
field with the cyclotron frequency v c — v+ + v-. Part (a) 
and (b) show the first and second half of the conversion, re- 
spectively. The solid circle in (a) indicates the initially pure 
magnetron motion. 
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FIG. 4: Excitation schemes: (a) standard excitation, (b) ex- 
citation with two 100 ms Ramsey pulses, (c) excitation with 
three 60 ms Ramsey pulses, and (d) excitation with four 40 
ms Ramsey pulses. The excitation amplitude is chosen in a 
way that at v c one full conversion from pure magnetron to 
pure cyclotron motion is obtained. Thus the sum of the grey 
colored areas is identical in all four schemes. 



B. Excitation schemes 

The precision in the determination of the cyclotron 
frequency strongly depends on the width of the cen- 
tral peak of the resonance. The narrower the reso- 
nance, i.e. the smaller the full-width-at-half-maximum 
(FWHM), the more precisely the center of the resonance 
can be determined. The pronounced central peak, which 
is clearly distinguishable from the usually smaller side 
bands, marks the point of the maximal conversion and 
thus of the maximal radial energy and the shortest time 
of flight from the trap to the detector. Ramsey's method 
of time-separated oscillatory fields promises to lead to 



narrower central peaks compared to the conventional ex- 
citation scheme. Therefore, a variety of different Ramsey 
excitation schemes were investigated with the aim of find- 
ing the one best suited for precision mass spectrometry. 
In [TlT | the line shapes are discussed in the quantum me- 
chanical formalism as probability distributions for the 
partial conversion of a given initial state into a state of cy- 
clotron motion. This general framework permits to take 
into account different assumptions on the initial state 
as well as statistical hypotheses, for example on phases. 
A more elementary approach consists in calculating the 
time development of the radius of the cyclotron motion 
R+(t) from the given initial conditions, using Eq. ©, ©, 
and ©. 
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1. The one-pulse excitation scheme 

In the conventional one-pulse excitation scheme the ex- 
ternal driving field is applied for a certain time interval 
t = Ttot with constant amplitude, as shown in Fig. |4] 
The lineshape represents the probability for the conver- 
sion of an initially pure state of magnetron motion into 
cyclotron motion as a function of the detuning 8 of the 
quadrupole field. For given values of r and g it is ob- 
tained in the general formalism of [ll| as 

F^r, 9) = ^sin 2 (^) . (9) 

Note that the maximum value 1 is reached when the 
driving field is at resonance, Vd = v Cl and the dura- 
tion of the excitation is chosen to satisfy the condition 
gr = (2n + 1)tt/2 with n = 0, 1, 2,.... In different words, 
at resonance (5 = 0) complete conversion of a pure mag- 
netron state into a state of pure cyclotron motion is 
achieved for excitation pulses of the duration of the 'con- 
version time' r c = ir/(2g) or odd multiples thereof. 

This result can be verified by application of Eq. ([5]) 
and ([7]) with the initial condition a+ (0) = and arbitrary 
a_(0), and computing the time development of the radius 
of modified cyclotron motion: 

Rl(r) = l^sin 2 (^).i? 2 _(0) = F 1 (S; r, 5 ).i? 2 _(0) . 

(10) 

With respect to time the shape of the excitation pulse is 
rectangular. Thus, in frequency space it is expected that 
the excitation resembles the intensity (i.e. the modulus 
squared) of the Fourier transform of a rectangular profile, 
namely (4g 2 /S 2 ) sin 2 (5r/2). The actual lineshape, how- 
ever, differs in two important respects from this Fourier 
transform: (a) At resonance (5 = 0) it is a periodic func- 
tion of r, describing the periodic conversion and recon- 
version of the magnetron and modified cyclotron modes, 
Fi(S = 0;t, g) = sin 2 (ttt / Y c ) , whereas the intensity of 
the Fourier transform increases proportionally to r 2 . (b) 
The central peak is actually narrower than for the Fourier 
transform of a rectangle. This can be deduced from the 
position 5q of the zero that separates the central peak 
from the first satellite peak, (Sot c ) 2 = 3ir 2 as compared 
to (^o T ) 2 = 47r 2 for the Fourier transform of the rectan- 
gle, assuming r = r c . 



2. The two-pulse excitation scheme and more general 
schemes 

Although it is not obvious on first sight, a close formal 
analogy exists between nuclear magnetic resonance on 
the one hand and the interconversion of the magnetron 
and cyclotron motional modes of an ion in a Penning trap 
due to quadrupole excitation on the other hand. This 
was shown in [1JJ using the concept of a Bloch vector. 



It is therefore reasonable to expect that the use of Ram- 
sey's method of separated oscillatory fields will lead to 
increased precision in mass spectrometry too. 

A symmetric n-pulse Ramsey cycle of a total duration 
r to t consists of n excitation intervals of duration n with 
(n— 1) waiting intervals of duration To in between, so that 
the total cycle time is T to t = nr\ + (n — 1)tq. Note that 
the n excitation pulses of the rf-quadrupole field must be 
coherent in phase. 

Assuming that the ion is initially in a state of pure 
magnetron motion, the probability for the conversion of 
magnetron quanta into quanta of the cyclotron motion 
by a 2-pulse Ramsey cycle with detuning 8 = tod — to c has 
been calculated to be 

F2(S;T ,Tx,g) = jcos (^Y^j sin (vrti) 

S . fSr { 
+ sm — 

UR V 2 / 

[C0S(CJ R T1)-1}} 2 . (11) 

This result and analogous ones for Ramsey excitation 
cycles with 3, 4, and 5 pulses can be found in [TT[ | . 
If the frequency of the quadrupole field equals the cy- 
clotron frequency v c and the amplitude of the field is 
chosen such that the n excitation intervals exactly add 
up to the conversion time t c , i.e. if the coupling con- 
stant g satisfies the relation m~i — t c = n/2g, then the 
profile function (flT|) reaches the value 1 at resonance, 
F 2 (6 = 0;T ,r c /2,g) = l. 

An elementary derivation of the 2-pulse profile function 
F 2 in Eq. dTTJ) is possible by applying Eq. ©, ©, and © 
three times successively to the time intervals < t < t\, 
ti < t < Ti + t , and ti + t < t < 2ti + t = T tot , with 
the initial condition ct + (0) = and arbitrary a_(0), in 
order to compute the endpoint of the time development 
of the radius of the cyclotron motion, i? + (T to t)- The first 
calculation yields the initial values a±(Ti) for the wait- 
ing period, the second calculation with g = yields the 
phase change during the waiting period and thus the ini- 
tial values ol±{ti + tq) for the second excitation period, 
the third step finally results in a±{2r\ + To) = a±(ji t) 
and thus in a result for the final radius of the cyclotron 
motion i£ + (T to t). Note that \d m Eq- © and © has 
to be replaced in the second and third step by the corre- 
sponding phases of the quadrupole field, c^ti + Xd and 
w d( T o + r i) + Xd, respectively. The final result is 

i?+(rtot) = F 2 (S;T ,n,g)-R 2 _(0). (12) 

The time development of the ion orbit during the two 
excitation periods is shown in Fig. [3^ and Fig.[5}D, during 
the waiting time the ion follows a rosette shaped orbit. 

Plotting the profile function (TrTj) with a fixed value of 
the waiting time To as a function of the detuning S one ob- 
tains the spectral lineshape of the 2-pulse Ramsey cycle. 
It bears a resemblance to the Fourier transform of a signal 
consisting of two rectangular pulses, but as for a single 
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FIG. 5: Percentage of converted quanta in the case of a 
quadrupolar excitation near ui c as a function of the frequency 
detuning 5' = S/2tt for different excitation schemes shown in 
Fig. [4] The total cycle time of all schemes is 300 ms. The line 
profile of a one-pulse excitation with 300 ms excitation time 
is shown in (a), (b) represents the line profile of a two-pulse 
excitation, each of 100 ms duration. The excitation time as 
well as the waiting time of the three-pulses excitation scheme 
in (c) is 60 ms. In (d) the four 45 ms pulses are interrupted 
by 40 ms waiting periods. The center frequency (5 — 0) is the 
cyclotron frequency u c for a given ion mass m. 



pulse there are important and characteristic differences. 
Lineshapes for higher order Ramsey cycles are calculated 
in an analogous fashion using the results of [HI]. In Fig. 
[5] generic results are displayed for n = 1, 2, 3, 4, assuming 
a total cycle time r to t = 300 ms and t\ = r c /n for all 
excitation schemes. Note that with the n-pulse excita- 
tion scheme a spectral distribution is obtained in which 
the valley between the first major sideband and the cen- 
tral peak contains (n — 2) small peaks, while the distance 
between the first major sideband and the central peak 
increases with n. 



III. IMPLEMENTATION AT THE MASS 
SPECTROMETER ISOLTRAP 

A. Experimental setup 



The triple-trap mass spectrometer ISOLTRAP 115 . 141 
installed at the on-line facility ISOLDE/CERN 0|l6| 
is dedicated to high-precision mass measurements of ra- 
dioactive nuclides. It reaches a relative mass uncertainty 
of below 10~ 8 [17]. A Paul trap and two Penning traps 
are the main parts of the apparatus as shown in Fig. [6] 
The ions from the continuous 60-keV radioactive beam of 
ISOLDE or from the stable alkali ion source are captured 
in-flight by a gas-filled linear Paul trap with segmented 
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FIG. 6: Schematic drawing of the mass spectrometer 
ISOLTRAP including the RFQ trap, the preparation and pre- 
cision Penning traps, as well as the reference ion source. In 
the inset a resonance for 39 Ca 19 F + ions with an excitation 
period of 900 ms is shown. To monitor the ion transfer and 
to record the time-of-flight resonance for the determination 
of the cyclotron frequency, micro-channel-plate (MCP) detec- 
tors or a channeltron are used. 



rods [18j]. Here, the ions are accumulated, cooled, and 
bunched. From the Paul trap the ions are guided, after 
passing two pulsed drift tubes in order to reduce the po- 
tential energy from 60 keV to about 100 eV, to a buffer- 
gas filled cylindrical Penning trap [l9j . This preparation 
trap is located in a superconducting magnet of 4.7-T field 
strength. Here, the ions are mass-selectively cooled [2(| 
with a resolving power of up to 10 5 . Thus, an isobar- 
ically clean ion cloud is obtained, which is transferred 
to the precision Penning trap. This hyperbolic Penning 
trap T4| |. which is placed in a second superconducting 
magnet of 5.9-T field strength with a field homogeneity of 
10~ 7 to 10~ 8 within one cubic centimeter and a temporal 
stability of ^ w 2 • 10~ 9 /h, provides the possibility to 
resolve even low-lying excited nuclear states [HI, H2] and 
serves for the actual mass measurement of nuclides with 
half-lives even below 100 ms (23l. l24j. 
In order to manipulate the ion motion in the precision 
trap by the application of external rf-fields, the ring elec- 
trode of the trap is four-fold segmented, as illustrated in 
Fig.O First, the magnetron radius of all ions is increased 
via a dipolar excitation at the magnetron frequency [25| . 
If required, there is the possibility to remove unwanted 
contaminations via a mass-selective dipolar excitation at 
the corresponding modified cyclotron frequency. Thereby 
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the ions' cyclotron radius is increased until the ions hit 
the electrode and are lost. A quadrupolar excitation at 
v c converts the magnetron motion into cyclotron motion. 
Finally, the ions are ejected out of the trap and pass the 
gradient of the magnetic field, which interacts with the 
orbital magnetic moment of the ions and accelerates them 
in the axial direction to the detector. This acceleration 
is proportional to the strength of the ions' magnetic mo- 
ment, i.e. to the radial energy obtained by the excita- 
tion. The time of flight after ejection from the trap to 
the detector is measured. This procedure is repeated for 
different frequencies of the quadrupolar excitation in the 
precision trap around the expected value of the cyclotron 
frequency. By the determination of the mean time of 
flight for the different excitation frequencies a time-of- 
fiight resonance curve is recorded (see inset of Fig. [S]). 
For appropriate excitation parameters the minimum time 
of flight is measured at the cyclotron frequency [2(| ■ For 
reference measurements, i.e. to calibrate the magnetic 
field, ions with well-known mass from a stable alkali ion 
source are used. 

The experimental standard deviation a{v c ) of the cy- 
clotron frequency v c is a function of the resolving power 
of the precision trap, i.e. the quadrupolar excitation 
time T q , and the total number N to t of recorded ions. 
The resolving power is Fourier limited by the duration 
of the quadrupolar excitation, which itself is limited by 
the half-life of the ion of interest in case of short-lived 
radionuclides. An empirical formula [13] describes this 
relation: 
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FIG. 7: Reduction of the width of the central peak due to 
the Ramsey method for constant cycle time rtot- For sym- 
metrical n-pulse excitation a Ramsey cycle is composed of n 
excitation pulses of duration ri = r c /n, where r c is the con- 
version time, separated by (n — 1) waiting periods of duration 
ro. Thus the total cycle time is r to t = t c + (n — 1)tq. The 
variable y = |<5Q n V to t/(7r\/3)| is a measure for the FWHM of 
the central peak for n-pulse excitation relative to the corre- 
sponding quantity for 1-pulse excitation, x — (n — l)ro/r to t 
is the fraction of the total cycle time spent during waiting 
periods. 



v c Vc \JNtot ■ T q 

where c is a dimensionless constant. In a large number 
of measurements with carbon clusters the constant c was 
determined for the ISOLTRAP mass spectrometer to be 
c = 0.898(8) 0. 



B. Reduction of the line-width 

Having discussed the standard quadrupolar excitation, 
what is the advantage of Ramsey's method of time- 
separated oscillatory fields? In many experimental sit- 
uations, the total time available for a complete mea- 
surement cycle has an upper limit, for example due to 
the lifetime of the radioactive species under investiga- 
tion. Thus, a precision gain simply by increasing the 
waiting time To may not be feasible due to experimen- 
tal limitations. Therefore, different excitation schemes 
are compared with respect to the predicted width of the 
central peak, assuming that a total time T to t is available 
to perform one complete Ramsey measuring cycle. For a 
symmetric n-pulse excitation it is r to t = nr\ + (n — 1)tq, 
where t\ denotes the duration of an excitation interval 
and To the duration of a waiting interval. For n = T c /n 
and an excitation at resonance a complete conversion of 



a pure magnetron state into a pure state of cyclotron mo- 
tion occurs, i.e. the highest possible degree of conversion. 
A suitable parameter to describe the width of the central 
peak is the position Sq 1 ^ of the zero separating the central 
peak from the first side band. For a one-pulse excitation 

(#o Ttot) 2 = 37r 2 . Therefore, for an n-pulse excitation 

the variable y = \6^T tot /(ny/3)\ = I^V^o^! represents 
the ratio of the zeros for n-pulse and 1-pulse excitation. 
To a very good approximation this variable equals the 
ratio of the full-width-half-maximum of the two excita- 
tion schemes. In Fig. [7] these ratios have been plotted 
for n = 2,3,4 as a function of x = (n — l)To/Ttot, i.e. 
the percentage of the total cycle time Ttot spent during 
waiting periods. From the graph it is obvious that with 
60% the two-pulse excitation scheme offers the largest 
width reduction relative to a one-pulse excitation. This 
excitation scheme is therefore favored for the application 
in high-precision mass spectrometry. 



C. Time-of- flight detection technique at 
ISOLTRAP 

As described above, the frequency determination via 
a time-of-flight detection technique [f| is based on the 
interaction of the magnetic moment of the orbiting ion 
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FIG. 8: The magnetic field amplitude from the precision trap 
to the end of the conversion section. Between the two marked 
points the amplitude decreases about 90% of its origin value. 



with the magnetic field gradient (Fig. [5]). Thus, for a 
detailed comparison with experimental data it is neces- 
sary to convert the theoretical line shapes for the energy 
conversion into time-of-flight spectra, which will be dis- 
cussed in the following. The kinetic energy in the radial 
motions is predominantly due to the cyclotron motion; 
the contribution of the magnetron mode is negligible be- 
cause of uj- <C 

pkin 



-m ■ (R\(T to t,Ud)u\ + R 2 _{T toU uj d )uj 2 _) 



-m ■ R + (T tat ,L0d)u + . 



(14) 



Here, R±{Ttot,^>d) are the radii of the two radial modes 
after the quadrupolar excitation has been applied with 
the frequency u>d- The magnetic moment of an ion with 
kinetic energy E r (u>d) m a magnetic field B = B ■ e z can 
be written as fl(w d ) = [E r (u)d)/B]e z . The interaction 
with the gradient of the magnetic field causes an axial 
force F z (u!d) — —fl(u)d) ■ V£? 2 on the ion, which leads 
to a reduction of the time of flight from the trap to the 
detector. This time of flight can be calculated with [26[ 



2 [Eq - q ■ V{z) - p(u d ) ■ B{z)) 



dz 



(15) 



where Eq is the total initial energy of the ion, V(z) 
and B(z) are the electric and magnetic fields, respec- 
tively, along the way from the trap to the detector. At 
Lo d = u c the magnetic moment is maximal and thus 
the time of flight minimal. Typical theoretical time- 
of-flight cyclotron-resonance curves for different excita- 
tion schemes using the radial energies calculated from 




FIG. 9: Calculated time-of-flight cyclotron resonances for dif- 
ferent excitation schemes. The total excitation and waiting 
time in the precision trap is 300 ms. The time of flight of 
a one-pulse excitation with 300 ms duration is shown in (a). 

(b) shows the time of flight of a two-pulse excitation, each 
of the pulses being 100 ms long. The three pulse excitation 

(c) is done by three 60-ms excitation periods and two 60-ms 
waiting periods. In (d) the four 45-ms pulses are interrupted 
by 40-ms waiting periods. 



the equations of conversion (see Eq. ([9]) and (| 1 1 1) and 
Ref. [ill]) are shown in Fig. [5] In each graph the time of 
flight of the ions from the trap to the detector is plotted 
as a function of the frequency detuning 5' = 5/2n with 
respect to the cyclotron frequency v c . 



IV. RESULTS 

To confirm the calculations, i.e. to determine the line- 
width reduction and, most importantly, to specify the 
precision gain due to the Ramsey excitation method, 
more than 300 time-of-flight resonance curves with differ- 
ent excitation schemes were recorded with the Penning 
trap mass spectrometer ISOLTRAP. The ion species for 
all measurements was 39 K + provided by the stable al- 
kali off-line ion source. Each resonance consists of about 
2500 ions in order to have identical statistics. To min- 
imize ion-ion interactions only time of flight measure- 
ments with at most five ions in the trap were taken into 
account. The efforts were concentrated on the specifi- 
cation of the uncertainty in the frequency determination 
for different excitation schemes. Figure [10] shows time- 
of-flight cyclotron resonances for the one-, two-, three-, 
and four-pulse excitation scheme. A fit of the theoreti- 
cally expected line shape (solid line) to the data points 
allows the determination of the FWHM and the cyclotron 
frequency v c along with its uncertainty 5v c . To perform 
these fits the standard evaluation program of ISOLTRAP 
[ItJ was extended in order to analyze the measured cy- 
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FIG. 10: The measured mean time of flight as a function of 
the quadrupolar excitation frequency with predicted curves 
fitted to the data. Here 39 K + ions from the stable alkali ion 
source were used, (a) Resonance of the conventional one-pulse 
excitation scheme with 300 ms excitation time, (b) a two-pulse 
excitation scheme with two times 100 ms excitation and 100 
ms waiting time, (c) three pulses, each 60 ms, interrupted by 
two waiting periods of 60 ms, (d) four 45 ms excitation pulses 
and three waiting periods of 40 ms. 



FIG. 11: The full-width-at-half-maximum (FWHM) values 
of the time-of-fiight cyclotron resonances as a function of the 
waiting time for the two-pulse excitation scheme is given. The 
total cycle times are r to t = 300 ms (squares), 600 ms (cir- 
cles), and 900 ms (triangles). The bold solid lines are the 
theoretically calculated values. Since the determination of 
the FWHM was performed manually from the fit curves, the 
error bars are conservatively estimated to be ±0.05 Hz. 



clotron resonances using the Ramsey method. 
The fit results concerning the FWHM are presented in 
Fig. [11] and Fig. [T^J Figure QT] shows results obtained 
with a two-pulse excitation scheme of different overall cy- 
cle times (r to t = 300 ms, 600 ms, 900 ms). The FWHM 
is given as a function of the waiting period. Due to field 
inhomogeneities and ion-ion interactions the data points 
are shifted slightly to higher FWHM values compared to 
theory. Similar results are shown in Fig. [T2] for differ- 
ent numbers of excitation pulses, where the total cycle 
time in the precision trap is constant r to t = 300 ms. The 
experimental values are on average 0.1 Hz higher than 
the theoretical ones. This line-broadening effect is due 
to the electric and magnetic field imperfections and ion- 
ion interactions, which were not taken into account in 
the calculations described above. A significant reduction 
of the FWHM for shorter excitation pulses with longer 
waiting periods in between can be observed. 
In Tab. [I] the experimental results are summarized, "cy- 
cle time" is the duration of the total cycle. The maxi- 
mum FWHM is the one of the standard one-pulse reso- 
nance curve. The minimum FWHM is measured using 
the longest possible waiting time, i.e. shortest possi- 
ble excitation time, which is determined by the maximal 
possible amplitude of the quadrupolar excitation field re- 
quired to obtain one full conversion from pure magnetron 
to pure cyclotron motion. The last column of Tab. Ogives 
the maximum line-width reduction (reduction gain= 
(max. FWHM - min. FWHM) /max. FWHM) using 
the different Ramsey excitation schemes. A remarkable 



TABLE I: The maximum and minimum experimental FWHM 
of the different excitation schemes for different cycle times 
are given. In addition the reduction gain is calculated. For 
further explanation see text. 



number 


cycle time 


max. FWHM 


min. FWHM 


reduction 


of pulses 


r to t 


Hz 


Hz 


gain % 


2 


300 


4.1 (0.1) 


2.6 (0.1) 


37.3 (2.0) 


2 


600 


2.1 (0.1) 


1.3 (0.1) 


38.1 (4.0) 


2 


900 


1.4 (0.1) 


0.9 (0.1) 


35.7 (5.9) 


3 


300 


4.1 (0.1) 


3.0 (0.1) 


26.5 (2.2) 


4 


300 


4.1 (0.1) 


3.3 (0.1) 


20.5 (2.7) 



reduction of close to 40% of the normal line-width is 
observed, similar to the results in the original work of 
Ramsey |8| . The reduction in line- width is especially im- 
portant in context with the achievable resolving power 
R = m/Am = v/Av. As theoretically predicted, the 
largest possible reduction is obtained by a two-pulse ex- 
citation scheme. However, the relative gain in reduction 
depends only weakly on the total cycle time r tot (see 
Tab. Hand Fig.Q3}. 

In Fig. [T3] the experimental uncertainty 8v c of the 
measured cyclotron frequencies for different numbers of 
pulses and different length of the total cycle time is plot- 
ted versus the waiting time t . Each data point repre- 
sents the mean value of three to ten individual measure- 
ments. The uncertainty 5v c decreases with increasing 
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FIG. 12: The full-width-half-maximum values of the time- 
of-flight cyclotron resonances as a function of the sum of all 
waiting periods To = r o f° r two, three and four excitation 
pulses. The total cycle time is r t ot=300 ms. The solid lines 
are the theoretically calculated FWHM values. 



waiting time. This is expected due to the decreasing 
FWHM at longer waiting times. If the cyclotron fre- 
quency uncertainty would only depend on the FWHM, 
a similar behavior as observed in Fig. [12] would be ex- 
pected. However, as mentioned before, it is obvious that 
there is also an effect of the overall line-shape, especially 
of the steepness of the curve and the pronounced side- 
bands, on the uncertainty in the frequency determination 
of v c . 

The excitation scheme used for the data points given in 
Fig. [T3] (a) consists of two, three, and four pulses, where 
the total excitation cycle is fixed to 300 ms. The un- 
certainty is obviously decreasing for shorter excitation 
pulses, i.e. longer waiting times To. In case of the two- 
pulse excitation it can be reduced by more than a factor 
of three (from 8v c « 27 mHz down to 8v c « 8 mHz) as 
compared to the conventional procedure just by chang- 
ing the excitation scheme to the two-pulse method. The 
uncertainty using the three- and four-pulse excitation 
scheme decreases to 8v c w 13 mHz and 5v c « 15 mHz, 
respectively. 

The result is summarized in Tab. HH where the max- 
imal and the minimal uncertainty for all investigated 
excitation schemes under identical experimental condi- 
tions are listed. As for the investigation of the FWHM, 
the scheme where the largest reduction in the cyclotron 
frequency uncertainty can be achieved is the two-pulse 
Ramsey scheme. Comparing the two-pulse excitation 
scheme with 300 ms, 600 ms, and 900 ms total cycle 
time, the tendency is similar to the results obtained for 
the FWHM of the resonances (see Fig. IT31 (b)). The 
slightly decreasing gain factor for two-pulse schemes with 
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FIG. 13: In (a) the uncertainty of the measured cyclotron 
frequency is given as a function of the waiting time for the 
two-pulse, three-pulse, and four-pulse excitation scheme. In 
(b) the uncertainty of the measured cyclotron frequency is 
given for the two-pulse scheme with a total cycle of 300 ms, 
600 ms, and 900 ms, where u — (n — l)ro/r to t is the fraction 
of the total cycle time spent during waiting periods. 



longer total cycle times must be assigned to the relative 
decrease of the statistical uncertainty in comparison to 
the constant systematic uncertainty. Unsymmetric exci- 
tation schemes where the individual Ramsey pulses have 
different lengths have also been investigated. However, 
the symmetric two-pulse excitation scheme remains the 
best one in respect to line-width reduction and uncer- 
tainty gain (see Tab. [TTJ) . 



V. FIRST ONLINE MASS SPECTROMETRY 
APPLICATION OF THE RAMSEY METHOD 

The first online mass measurement by the Ramsey 
excitation method was carried out for the short-lived 
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TABLE II: The maximum and the minimum experimental 
uncertainties of different excitation schemes and excitation 
times are listed. The improvement factor is given in the last 
column. 



number 


cycle time 


max. 


min. 


Improvement 


of pulses 


Hot 


uncertainty 


uncertainty 


factor 






Hz 


Hz 




2 


300 


0.027 


0.008 


3.4 (0.4) 


2 


600 


0.017 


0.006 


2.8 (0.3) 


2 


900 


0.010 


0.004 


2.5 (0.3) 


3 


300 


0.027 


0.013 


2.1 (0.2) 


4 


300 


0.027 


0.015 


1.8 (0.2) 
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FIG. 14: Time of flight for 39 Ca 19 F+. A two-pulse Ramsey 
scheme was chosen with two 100 ms duration excitation peri- 
ods interrupted by a 700 ms waiting period. The solid curve 
is a fit of the theoretically expected line shape to the data. 

nuclides 38 Ca and 39 Ca, which have half- lives of only 
T 1/2 ( 38 Ca) = 440 ms and T 1/2 ( 39 Ca) = 860 ms, respec- 
tively [23. In order to suppress 38 K + contaminations, 
the 38 Ca ions were delivered in form of the sideband 
molecule 38 Ca 19 F+. In the inset of Fig. [6] a resonance 
curve of 39 Ca 19 F + is shown for which the ions were ex- 
posed to a continous quadrupolar radiofrequency excita- 
tion of 900 ms duration. The cyclotron frequency of this 
resonance has been determined with an uncertainty of 
5v c = 32 mHz. Fig. [14] shows a resonance of the same 
species with the same number of collected ions 2500), 



where a two-pulse Ramsey excitation scheme was used. 
The two excitation pulses had a duration of 100 ms in- 
terrupted by a waiting period of 700 ms. Thus, the total 
time Ttot = 900 ms for which the ions remained in the 
trap was identical. Here, the cyclotron frequency was 
determined with an uncertainty of only 5v c = 9 mHz. 
In comparison, the statistical error in the frequency de- 
termination could be reduced by more than a factor of 
three, keeping the number of ions and the storage time 
constant. This is a tremendous gain factor, especially for 
mass measurements on short-lived radionuclides since the 
required measurement time to reach a certain statistical 
uncertainty can be reduced by almost an order of mag- 
nitude. 



VI. CONCLUSIONS AND OUTLOOK 

The experimental studies described in this paper 
demonstrate that Ramsey's method of time-separated os- 
cillatory fields can be applied to excite the ion motion in 
a Penning trap. The Ramsey technique improves signif- 
icantly the statistical uncertainty in high-precision mass 
spectrometry on short-lived radionuclides. We performed 
systematic experimental investigations for different ex- 
citation patterns with two, three, and four excitation 
pulses. We observed the expected reduction of the line 
width of almost a factor of two along with a gain in pre- 
cision in the frequency determination. This leads to an 
important gain in precision of the frequency determina- 
tion. The new findings were demonstrated for a stable 
nuclide as well as for a short-lived radionuclide in an on- 
line measurement. 

An optimized Ramsey excitation scheme with two 
pulses of short duration interrupted by a long waiting 
period results in an improvement of the statistical un- 
certainty in the cyclotron-frequency determination by 
more than a factor of three compared to the conventional 
scheme, without any further experimental changes as e.g. 
the number of detected ions for a resonance or reduction 
of scan detuning width. Since the Ramsey method opens 
a door to higher precision in Penning trap mass spec- 
trometry, the application in other Penning trap setups 
is already in preparation and under investigation, as e.g. 
at SMILETRAP (Stockholm) using highly-charged sta- 
ble ions [Io| or at SHIPTRAP (GSI, Darmstadt) using 
short-lived fission fragments ^29] . It can be expected that 
the Ramsey method will find a wide-spread application 
in high-precision mass spectrometry of atomic nuclei. 
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